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EXECUTIVE  SUMMARY 


New  Large  Aircraft  (NLA)  pose  unique  firefighting  challenges  not  previously  considered  for 
traditional  aircraft.  Specifically,  questions  have  arisen  regarding  the  applicability  of  current 
firefighting  protection  standards  to  nonconventional  design  changes,  such  as  fuselage  shape 
modifications,  enhanced  material  compositions,  new  fuel  storage  locations,  and  unique  passenger 
loading  configurations.  To  address  fuselage  shape  concerns,  a  computational  fluid  dynamic 
(CFD)  modeling  strategy  was  developed  to  quantify  heat  transfer  to  the  aircraft  for  a  given 
aircraft  geometry  and  hydrocarbon  pool  fire  condition. 

A  brief  review  on  large-scale  pool  fires  and  their  interactions  with  engulfed  bodies  both  from  an 
experimental  as  well  as  computational  perspective  is  initially  presented.  Experimental  data 
published  by  the  Fire  Science  and  Technology  Division  of  Sandia  National  Laboratories  and  the 
Naval  Air  Warfare  Center  Weapons  Division  on  a  large-scale  mock  fuselage  engulfed  in  a 
hydrocarbon  pool  was  used  to  validate  the  CFD  methodology  for  low-  to  medium  crosswind 
conditions;  a  range  defined  by  combustion  and  turbulence  modeling  assumptions.  The  validated 
CFD  process  was  then  applied  to  two  NLA  geometries  (Boeing  777  and  Airbus  A380)  and  one 
traditional  platform  (Boeing  707)  to  exercise  the  model’s  applicability  using  the  National  Fire 
Protection  Association  403  standard  practical  critical  area  fire  condition. 

The  CFD  predictions  supported  experimental  conclusions  showing  increased  crosswinds  amplify 
heat  transfer  to  the  aircraft  surface  due  to  enhanced  turbulent  fuel-air  mixing.  Based  upon  the 
predicted  CFD  fire  plume  structure  and  aircraft  surface  heat  transfer  magnitudes,  nominal 
changes  in  aircraft  geometry  exposed  to  similar  scale  flame  and  atmospheric  conditions  pose  no 
extraordinary  firefighting  challenge. 
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(a)  (b) 

Figure  2.  The  NLA  Mockup  (a)  Engulfed  in  a  JP-8  Pool  Fire  and  (b)  After  the  JP-8  Pool  Fire 

Using  Infrared  Imagery 


Although  full-scale  fire  testing  is  critical  to  fire  research,  it  has  some  disadvantages.  Fuel  costs 
can  be  prohibitive.  Outdoor  studies,  which  are  necessary  for  NLA  analysis  due  to  the  massive 
spatial  scales  involved,  must  contend  with  chaotic  atmospheric  influences  that  typically  result  in 
significant  variance  between  repeated  trials.  Experimental  setups  also  encounter  data  acquisition 
(DAQ)  challenges  regarding  measurement  certainty;  and  embedded  instrumentation  is  often 
restricted  to  a  few  finite  points,  offering  a  limited  scope  of  the  flow  field.  In  addition,  the 
destructive  combustion  environment  results  in  reduced  life  expectancy  for  both  fabricated  test 
and  DAQ  equipment. 

1.2  SCOPE. 

To  supplement  full-scale  fire  tests,  a  computational  fluid  dynamic  (CFD)  simulation  strategy  was 
developed.  The  objective  of  the  CFD  was  to  predict  the  major  combustion  plume  flow  features 
as  well  as  thermal  heat  transfer  magnitudes  to  aircraft  bodies  immersed  in  hydrocarbon  pool  fires 
that  are  consistent  with  an  accidental  post-crash  fire  event.  The  CFD  offers  an  alternative 
approach  for  obtaining  predictive  estimates  of  fire  scenarios  too  impractical  or  costly  to  evaluate 
experimentally. 

This  research  consisted  of  a  literature  review,  selection  of  the  CFD  software  suite,  and  simulation 
development.  The  literature  survey  briefly  introduces  the  physical  principles  associated  with 
large-scale  hydrocarbon  liquid  pool  fires,  their  interaction  with  engulfed  bodies,  along  with  a 
brief  overview  of  pool  fire  modeling. 

An  industry- accepted,  commercially  available  CFD  software  suite  entitled  ANSYS  FLUENT® 
was  selected  as  the  analysis  tool  to  ensure  continued  physical  model  development  and 
accessibility  in  the  public  domain.  Its  hybrid  mesh  capability  also  provides  advantages  involving 
complex  geometries  in  which  most  exclusive  fire  modeling  codes  still  face  limitations.  The 
methodology  is  initially  presented  and  validated  using  experimental  data  published  by  Sandia 
National  Laboratories  (SNL)  in  conjunction  with  the  Naval  Air  Warfare  Center  Weapons 
Division  (NAWCWPNS),  and  then  applied  to  particular  aircraft  of  interest. 
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The  primary  objective  of  this  research  was  to  develop  the  CFD  simulation.  Special  attention  was 
paid  to  the  role  aircraft  geometry  plays  in  altering  localized  pool  fire  combustion  flow  dynamics, 
specifically  fuselage  shape,  wing  planform,  and  engine  configuration.  Using  similar  fire 
conditions,  results  were  obtained  comparing  the  time-averaged  thermal  response  between 
different  aircraft.  Results  were  intended  to  provide  supplemental  risk  assessment  metrics  to 
support  firefighting  protection  needs  for  current  and  future  generation  aircraft.  The  long-term 
contribution  of  this  work  aims  to  support  an  integrated  modeling  simulation  strategy  accounting 
for  aircraft  post-crash  dynamics,  accidental  fire  scenarios,  as  well  as  fire  suppression  interaction 
to  supplement  auxiliary  research  initiatives. 

2,  LITERATURE  REVIEW. 

2.1  OVERVIEW. 

Literature  was  compiled  that  examined  the  flow  characteristics  of  large-scale  hydrocarbon  pool 
fires  described  experimentally  as  well  as  computationally,  especially  those  involving  interactions 
with  immersed  bodies.  For  this  study,  a  large-scale  pool  fire  is  defined  as  a  non-premixed, 
buoyant  diffusion  flame  established  over  a  two-dimensional  (2-D)  horizontal  liquid,  heavy 
hydrocarbon  fuel  surface.  A  heavy  hydrocarbon  is  characterized  by  properties  consistent  with 
aviation  jet  fuels  representing  common  diesel  derivatives,  such  as  Jet  A  and  JP-8.  Pool 
diameters  are  assumed  to  be  much  greater  than  1  m  to  coincide  with  spatial  scales  associated 
with  full-size  commercial  aircraft. 

Due  to  the  lack  of  computational  power  available  to  solve  all-inclusive  field  models  with 
accuracy  and  resolution,  the  majority  of  relevant  information  prior  to  the  mid-1990s  focused  on 
experimental  and  observational  data.  The  following  summary  focuses  on  three  categories 
including  the  physical  characteristics  of  large-scale  pool  fire  dynamics,  their  interaction  with 
engulfed  bodies,  and  the  computational  developments  taken  to  complement  analytical  and 
experimental  efforts. 

2.2  PHYSICAL  CHARACTERISTICS  OF  LARGE-SCALE  POOL  FIRES. 

Several  reviews  discussing  hydrocarbon  pool  fire  physics  summarize  most  of  the  literature 
compiled  over  the  past  few  decades  [6-10],  most  recently  by  Steinhaus  et  al.  [11],  who  focuses 
on  large-scale  pool  fires  resulting  from  accidental  spills.  Steinhaus  separates  the  fire  into  five 
distinct  zones  including  (1)  the  liquid  fuel,  (2)  the  unbumed  fuel  vapor  cone  immediately  above 
the  liquid  fuel,  (3)  a  luminous  combustion  flame  region  surrounding  the  fuel  vapor  cone,  (4)  an 
intermittent  turbulent  combustion  zone  above  the  luminous  flame,  and  (5)  a  nonreacting  buoyant 
plume  shrouding  the  interior  zones  [11],  These  zones  are  described  as  being  physically 
interdependent  by  nature  requiring  mathematical  closure  to  provide  quantitative  definition  to 
each.  Together,  the  system  harbors  complexities  associated  with  a  multiscale,  multiphase, 
chemically  reacting  turbulent  environment  [12],  Steinhaus  recaps  metrics  critical  to 
characterizing  the  flow  through  each  zone:  fuel  bum  rate  (or  mass  loss  rate),  heat  release  rate  (or 
the  total  energy  released  from  the  fire),  flame  height  (loosely  defined  as  distance  from  the  fire 
base  to  the  flame  tip  with  at  least  50%  intermittency),  smoke  production  rate,  and  radiation 
released  (or  emissive  power).  These  factors  are  driven  by  a  variety  of  elements  ranging  in 
sophistication.  A  few  common  descriptors  include  pool  geometry,  fuel  composition,  ventilation 
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conditions,  adjacent  geometry,  and  bounding  substrate  [11].  Experimental  methods  quantifying 
these  characteristics  are  uniquely  discussed  by  several  researchers  employing  various  techniques 
[13-23], 

It  is  speculated  that  the  principal  vortical  flow  structures  that  uniquely  identify  large-scale  pool 
fires  similar  to  that  shown  in  figure  2  are  the  result  of  baroclinic,  or  stratification  of  vortex  pairs 
produced  from  the  buoyant  collision  of  sequentially  smaller  flow  structures  [12].  Generated  by 
the  turbulent  energy  cascade,  pool  fires  contain  a  wide  range  of  temporal  and  spatial  scales 
within  these  structures  responsible  for  various  characteristics.  Larger  turbulent  scales  at 
relatively  low  frequencies  (often  described  as  an  intermittent  pulsation,  or  vortical  eddy 
shedding)  shape  overall  plume  dynamics,  which  drives  global  air  ingestion  levels  and  low-order 
interaction  between  the  combustion  plume  and  surrounding  air.  Conversely,  finer,  high- 
frequency  turbulent  flow  scales  directly  control  local  fuel-air  mixing  and  combustion  rates  [10]. 
Visible  pulsation  frequencies  dominated  by  large-scale  structures  are  well  described  and  are 
found  to  decrease  with  pool  fire  diameter  through  simple  empirical  correlations  using  Froude  and 
Strouhal  dimensionless  numbers.  These  transient  pulsations  are  also  believed  to  be  highly  reliant 
upon  soot  production  [12]. 

Large-scale  pool  fire  heat  transfer  is  dominated  by  radiation  where  the  fuel  bum  or  recession  rate 
becomes  nearly  independent  of  increasing  pool  diameter.  This  is  due  to  saturated  soot  levels 
creating  an  optically  thick  environment  resulting  in  strong  radiation  feedback  to  the  liquid  fuel- 
vapor  surface.  In  opposition,  small-scale  (<2  m)  pool  fires  are  dominated  by  convection  heat 
transfer  where  soot  production  is  sustained  at  lower  levels.  Ventilation  or  wind  effects  can  also 
greatly  influence  large-scale  dynamics  through  enhanced  mixing  in  the  combustion  region, 
which  causes  increased  chemical  reaction  rates  and  elevated  flame  temperatures.  In  addition, 
high  crosswind  environments  skew  radiation  profiles  that  can  significantly  change  fuel 
vaporization  levels,  which  further  complicates  the  heat  transfer  process.  Pool  bounding 
structures  can  also  alter  pool  fire  dynamics  by  modifying  thermal  resistance.  Soot  production 
remains  a  challenging  but  important  subject  matter  as  its  production  and  transport  is  heavily 
reliant  upon  high-order,  nonlinear  interacting  turbulence,  temperature,  and  fuel  composition 
effects  [12].  Although  various  studies  have  characterized  different  physical  aspects  of  soot,  it 
remains  a  focal  point  of  ongoing  research  [24-26]. 

2.3  OBJECTS  IMMERSED  IN  LARGE-SCALE  POOL  FIRES. 

Pool  fire  physics  are  further  complicated  by  the  introduction  of  engulfed  objects  to  the  flame 
front.  These  object-thermal  interactions  increase  in  significance  as  the  ratio  of  object  size  to  fire 
size  increases.  The  majority  of  experimental  work  conducted  on  objects  engulfed  in  large-scale 
hydrocarbon  pool  fires  stems  from  joint  investigations  lead  by  SNL  and  NAWCWPNS. 
Research  sponsored  by  the  U.S.  Department  of  Energy  (DOE),  the  Nuclear  Regulatory 
Commission,  and  the  Department  of  Transportation  examined  the  effect  hydrocarbon  pool  fire 
dynamics  on  engulfed  mockup  structures  simulating  the  aircraft  post-crash  fire  environment. 
Similar  studies  pertained  to  shipping  containers  of  hazardous  materials  in  similar  combustion 
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environments  [27  and  28].  Objectives  from  these  studies  summarized  by  Koski  et  al.  are  listed 
below. 


Enhance  the  understanding  of  fire  phenomenology. 

Estimate  empirical  parameters  for  Risk  Assessment  Compatible  Fire  Models  (RACFMs). 
Provide  validation  assistance  for  fire  models. 

Expand  the  thermal  response  database  of  objects  engulfed  in  pool  fires  [29], 

Early  experiments  conducted  on  a  water-cooled  (constant  surface  temperature)  calorimeter 
helped  determine  radiation  absorption  coefficients  and  demonstrated  that  gray  gas  models 
provide  a  reasonable  description  in  a  soot-laden  environment  [29],  The  first  in  a  series  of  mock 
fuselage  tests  involving  large  steel  cylinders  engulfed  in  jet  fuel  fires  quantified  the  effects  of 
turbulent  combustion  environments  exposed  to  various  conditions.  Results  showed  a  rise  in 
flame  temperature  and  surface  heat  flux  accompany  turbulent  recirculation  zones  subject  to 
accelerated  fuel-air  mixing.  These  elevated  zones  were  found  on  the  leeward  edge,  or  wake 
region  of  the  immersed  object  with  gains  proportional  to  the  crosswind  strength  [30]. 
Subsequent  SNL  testing  analyzing  object-thermal  response  relative  to  pool  fire  diameter  can  be 
found  in  associated  reports  [31-33].  These  reports  focus  on  instrumentation  development  and 
error  quantification.  Data  from  the  test  series  published  by  Gritzo  and  Suo-Antilla  [33]  was 
extracted  to  provide  validation  data  for  the  CFD  modeling  approach  employed  for  this  study. 

2.4  LARGE-SCALE  POOL  FIRE  MODELING. 

The  majority  of  relevant  experimental  work  supports  the  verification  and  validation  of 
computational  fire  models  to  aid  in  probabilistic  risk  assessment.  In  the  mid-1990s,  the  FAA 
commissioned  SNL  via  an  interagency  agreement  to  analyze  the  technical  issues  involved  with 
fire  modeling,  as  well  as  to  plan  research  steps  to  address  future  requirements.  Fire  models  were 
broken  down  into  four  major  categories  [5], 

Zone  models  that  are  largely  based  upon  correlated  experimental  data  that  was  tuned  for 
specific  regimes  typically  providing  rapid  solutions  for  compartmental  or  enclosure  fires. 

Field  models  (the  application  of  this  effort)  focused  on  first  principle,  physical 
interactions  that  generally  require  intensive  computing  power,  are  more  robust  than  zonal 
models,  and  are  less  restrictive  on  geometry  driven  scenarios  of  interest. 

Pyrolysis,  or  material  burning  models,  that  require  closed-loop  feedback  between  gaseous 
combustion  and  solid  devolatilization/breakdown. 

RACFM  derived  from  a  combination  of  the  above  models  that  was  optimized  to  solve  the 
dominant  physical  processes  driving  combustion  and  heat  transfer  [5]. 


SNL  concluded  that  field  models  offered  the  most  reliable  approach  for  fire  modeling,  noting 
further  work  in  turbulence  modeling,  soot  generation,  and  radiation  properties  of  soot-laden 
gases  was  needed.  Model  sensitivity  analyses  and  experimental  validation  data  were  highlighted 
as  key  contributors  to  an  overall  proposed  research  plan.  Two  SNL  fire  research  codes  analyzed, 
VULCAN  and  the  SNL  Gray  Gas  Model,  provided  agreeable  results  compared  to  available 
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experimental  data.  In  general,  it  was  found  that  computational  fire-modeling  tools  provided 
valuable  qualitative  and  quantitative  analysis  of  aircraft- fire-related  problems  [5], 

Because  hydrocarbon  pool  fire  dynamics  are  dominated  by  a  wide  range  of  length  and  time 
scales,  high-fidelity  resolution  of  these  interacting  scales  pose  the  greatest  challenge.  Because 
they  are  computationally  efficient,  physical  modeling  approaches  resolve  the  dominant  physical 
transport  process  typically  defined  by  a  finite,  narrow  range  of  scales.  The  wider  the  spectrum, 
the  more  computationally  intensive  the  solution  process  becomes.  A  review  of  all  the  interacting 
fluid  and  combusting  temporal  and  spatial  scales  can  be  found  in  a  summary  by  Tieszen  [34],  A 
detailed  overview  of  several  industry-accepted  fire  field  submodels  and  their  advantages  and 
disadvantages  can  be  found  in  a  text  by  Yeoh  and  Yuen  [35].  A  short  survey  of  supplemental 
work  investigating  various  numerical  approaches  to  solving  different  aspects  of  the  general 
hydrocarbon  pool  fire  problem  is  also  highlighted  [36-40].  The  physical  submodels 
implemented  for  this  study  are  briefly  discussed  in  section  3.3. 

Fire  Dynamics  Simulator  (FDS),  a  popular  field-modeling  approach  available  through  the 
National  Institute  of  Standards  and  Technology  (NIST),  is  a  CFD  solver  specifically  formulated 
for  low-speed,  thermally  driven  flows  characteristic  of  accidental  fires.  Although  it  is  considered 
a  research  code,  several  add-on  packages  tailored  for  flow  visualization  and  user-friendly  post¬ 
processing  make  FDS  an  attractive  computational  engineering  tool  for  fire  modeling.  For  the 
most  part,  the  developmental  history  of  FDS  focused  on  internal  fire  scenarios;  however 
continued  development  has  extended  applicability  to  external  hydrocarbon-fueled  environments. 
Mathematical  model  details  and  user  implementation  can  be  found  at  the  official  NIST  website 
[41].  FDS  currently  relies  strictly  on  rectilinear  coordinates,  or  the  enforcement  of  structured 
grid  topologies  in  the  geometric  domain.  Because  high-fidelity  resolution  of  the  complex 
curvature  of  commercial  aircraft  geometries  required  unstructured  mesh  entities  for  this 
particular  study,  FDS  was  not  considered. 

More  recent  research  by  the  Center  for  the  Simulation  of  Accidental  Fires  and  Explosions 
(CSAFE),  a  DOE  Advanced  Simulation  and  Computing  Center,  has  further  advanced 
understanding  of  large-scale  pool  fire  dynamics  through  the  development  of  a  unique, 
sophisticated  in-house  simulation  framework  in  conjunction  with  dedicated  computational 
resources.  Their  work,  similar  to  SNL,  centers  on  risk  assessment  of  objects  immersed  in  pool 
fires  from  the  perspective  of  weapons  safety.  Although  the  program  was  decommissioned  in 
November  2010,  uncertainty  quantification,  radiation  modeling  in  complex  fire  simulations,  and 
advanced,  detailed  chemistry  methods  for  soot  generation  and  deposition  are  among  the  many 
critical  subfields  still  actively  investigated.  A  complete  list  of  associated  CSAFE  publications 
can  be  found  at  the  official  CSAFE  website  [42],  in  addition  to  synopsis  papers  summarizing 
past  work  [43-45], 

3.  COMPUTATIONAL  MODEL  VALID ATION. 

3.1  OVERVIEW. 

Validation  case  studies  were  selected  based  upon  well-established  test  procedures,  clarity  of 
overall  design  and  boundedness,  as  well  as  depth  of  measurement  detail.  Geometric  simplicity 
was  also  a  contributing  factor  to  reduce  errors  prone  to  complex  meshing  environments.  CFD 
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models  were  constructed  to  reflect  the  experimental  spatial  domain  and  appropriate  boundary 
conditions  to  compare  results.  A  test  sequence  reported  by  the  Fire  Science  and  Technology 
Division  SNL  in  coordination  with  the  NAWCWPNS  was  selected.  The  test  series  was 
established  to  support  Weapons  System  Safety  Assessments  in  the  scenario  of  transport  aircraft 
exposed  to  a  large  fuel  fire.  Its  main  focus  was  the  study  of  heat  transfer  effects  of  a  cylindrical 
fuselage  shape  reminiscent  of  a  military  transport  aircraft  partially  engulfed  in  a  hydrocarbon 
pool  fire.  Test  objectives  mandated  the  data  be  presented  in  a  context  suitable  for  comparison  to 
computational  model  predictions.  Its  purpose  was  intended  to  spur  further  development  of  fire 
models  that  could  be  used  to  address  other  complex  fire  scenarios  of  interest.  In  addition,  the 
study  was  part  of  a  minority  of  experiments  that  detail  thermal  response  at  similar  physical  scales 
to  that  of  commercial  passenger  aircraft  where  engulfed  body  and  flame  dynamics  are  strongly 
interdependent  [33], 

3.2  EXPERIMENTAL  SUMMARY. 

Computational  validation  models  were  defined  to  reflect  setup  and  operating  conditions  from 
experimental  trials  conducted  at  the  NAWCWPNS  CT-4  test  site.  Test  geometry  consisted  of  a 
corrugated  3.66-m  diameter,  18.28m-long  culvert  welded  shut  at  both  ends  and  supported  0.6  m 
above  the  leeward  edge  of  an  18.9-m  diameter  JP-8  jet  fuel  pit.  The  test  pit  was  25-cm  deep, 
about  10  cm  of  JP-8  floating  atop  approximately  15  cm  of  water.  Ignition  was  prompted  at  3 
locations  via  a  110  V  electrical  signal.  Testing  concluded  once  all  the  fuel  was  burned.  The 
cylindrical  culvert  was  made  of  16-gage  (0.16-cm-thick)  carbon  steel  wrapped  internally  with 
approximately  5  cm  of  ceramic  fiber  blanket  insulation  to  protect  instrumentation  and  prolong 
structural  life.  An  array  of  heat  flux  gages  (HFGs)  and  thermocouples  (TCs)  measuring  total 
heat  flux  and  temperature,  respectively,  were  arranged  around  the  surface  of  the  cylinder, 
beneath  and  along  the  liquid  fuel  surface,  and  in  suspended  vertical  locations  above  the  liquid 
fuel  surface.  The  cylinder  surface  was  characterized  in  three  locations  equally  spaced  3.66  m 
about  the  culvert  centerline  identified  as  the  central  (C),  western  (W),  and  eastern  (E)  location. 
Similar  instrumentation  was  also  placed  just  outside  of  the  flame  zone  to  compare  relative 
effects.  Differential  pressure  gages  (DPGs)  were  installed  lengthwise  on  the  culvert  surface  on 
the  windward  and  leeward  sides  parallel  to  the  prevailing  wind  direction  at  the  central,  western, 
and  eastern  locations.  A  final  set  of  DPGs  was  stationed  on  the  fuel  pool  surface  and  outside  of 
the  flame  zone  to  monitor  pressure  gradients  generated  by  the  combustion  process  relative  to  the 
far  field.  Free  stream,  time-accurate  wind  speed  and  direction  were  recorded  using  vane-type 
gages  positioned  far  enough  away  from  the  fire  plume  to  minimize  exposure  to  air  entrainment 
induced  by  the  fire  plume  [33].  Figure  3  depicts  the  actual  experimental  setup  with  an  aerial 
schematic  view  highlighting  the  approximate  locations  of  HFG,  TC,  and  DPG. 
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Figure  3.  The  Experimental  Set-up  Showing  (a)  the  C-141  Mock  Fuselage  in  a  20-m-Diameter 
Fuel  Pit  and  (b)  an  Overhead  Schematic  With  Approximate  DPG,  HFG,  and  TC  Locations 

(Adapted  from  reference  33.) 


During  the  DAQ  process,  periods  of  quasi-steady  state  crosswind  and  fire  plume  behavior  were 
recorded.  Three  experiments  out  of  a  possible  eight  reported  were  identified  for  comparison  to 
CFD  model  results  covering  what  will  be  referred  to  as  a  low  (1.9  m-s'1),  medium  (5.4  m-s'1), 
and  high  (10.2  m-s'1)  mean  velocity  crosswind  condition.  These  cases  were  selected  because 
they  had  the  lowest  root  mean  square  (RMS)  velocity  magnitude  fluctuations,  or  most  consistent 
wind  speed  and  direction  out  of  the  available  trials.  Table  1  summarizes  the  three  experimental 
cases  selected  for  comparison  to  CFD  results.  Time  after  ignition  denotes  the  quasi-steady  time 
period  where  DAQ  took  place.  Wind  velocity  and  direction  are  described  by  an  average 
magnitude  and  its  associated  standard  deviation. 


Table  1.  A  Summary  of  Experimental  Case  Studies  Employed  for  Computational 

Model  Validation  [33] 


CaseName 

Time  Afterignition 
(s) 

Wind 

V  elocityMagnitude 
(rms'1) 

Wind  Direction 

(°) 

Low 

225-350 

1.9  ±0.2 

-36.9  ±5.7° 

Medium 

400-575 

5.4  ±1.2 

11.4  ±12.5° 

High 

300-600 

10.2  ±1.7 

-22.7  ±8.3° 

An  average  fuel  recession  rate  of  4.4  mmxnirf1  within  the  fuel  pit  was  derived  from  the  unsteady 
temperature  change  recorded  by  TCs  located  along  vertical  arrays  positioned  through  the  water, 
fuel,  and  flame  zone.  A  fuel  consumption  rate  of  0.058  kgxn'2^'1  was  extrapolated  based  upon 
fuel  pool  geometry  and  physical  properties  associated  with  JP-8.  An  average  fuel  surface 
temperature  of  508  K  was  recorded  for  all  trials.  Figure  4  shows  instantaneous  screen  captures 
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of  a  fully  developed  fuel  fire  engulfing  the  mock  fuselage  at  the  medium  crosswind  condition 
and  the  medium  crosswind  time-averaged  over  the  quasi-steady  period  [33], 


(a)  (b) 

Figure  4.  Adapted  Instantaneous  Screen  Capture  of  the  Mock  Fuselage  (a)  Fully  Engulfed  in 
Fire  Exposed  to  a  Medium  Wind  Condition  and  (b)  a  Time-Averaged  Video  of  the  Same 

Conditions  (Adapted  from  reference  33.) 

3.3  COMPUTATIONAL  MODEL  SETUP. 

It  is  well  established  that  multiscale  turbulent  flow  mixing,  interacting  chemical  kinetic 
combustion  processes,  radiation  heat  transfer  effects,  and  soot  production  are  the  main  features 
that  drive  large-scale  accidental  fire  plume  dynamics  [1 1  and  12].  Multiple  physical  submodels 
were  solved  simultaneously  using  ANSYS®  FLUENT  to  approximate  this  complex 
phenomenology.  Through  trial  and  error  of  mechanisms  with  varying  degrees  of  fidelity,  a 
solution  framework  was  created  to  accurately  and  practically  predict  the  experimental  validation 
conditions.  The  software’s  compatibility  within  a  multiphysics  environment  and  support  for 
complex  mesh  topologies  reinforced  its  choice  over  alternatives.  Table  2  summarizes  the 
computational  methodology  implemented  with  each  physical  submodel,  which  are  summarized 
in  subsequent  sections.  The  approach  is  applied  to  simulate  the  fire  environments  measured  in 
the  experimental  validation  study  presented  in  section  3.2.  Geometric  mesh  and  boundary 
condition  details  conclude  the  summary. 

Table  2.  A  Summary  of  the  Computational  Modeling  Approach  Taken  to  Predict  Aircraft  Crash 

Pool  Fire  Events 


Physical  Submodel 

Physical  Submodel  Description 

3-D  Navier  Stokes 

Conservation  equations  of  fluid  motion 

Non-premixed  combustion 

Chemical  reaction  kinetics 

Realizable  k- 

Turbulence 

Discrete  ordinates  (DO) 

Radiation  heat  transfer 

Moss-Brookes 

Soot  production 
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3.3.1  Turbulence. 


Large-diameter  (>1  m)  hydrocarbon  pool  fires  are  characterized  by  their  periodically  oscillating 
but  intrinsically  chaotic  thermally  buoyant  plume.  These  physical  traits  are  indicative  of  a  fully 
turbulent  flow  structure.  Turbulent  flow  quality  is  defined  by  the  Reynolds  number  (Re),  which 
describes  the  ratio  of  inertial  to  viscous  forces  in  a  particular  system;  the  higher  the  Re,  the  more 
turbulent  the  flow  quality.  Another  metric  examines  the  turbulence  intensity  (TI)  contained 
within  the  flow,  or  ratio  of  magnitudes  between  the  fluctuating  and  mean  flow  components. 
Table  3  outlines  the  TI  and  the  Re  associated  with  each  validation  case.  Because  Re  designation 
relies  on  a  characteristic  length,  two  values  were  calculated — one  based  upon  the  fuel  pit 
diameter  (18.9  m),  and  another  based  upon  the  mock  fuselage  diameter  (3.66  m).  Transition  to 
turbulent  flow  over  a  2-D  cylinder  (i.e.,  mock  fuselage)  begins  around  a  Re  of  105  [46].  With  the 
lowest  calculated  value  from  table  3  on  the  order  of  1 05,  the  flow  domain  was  assumed  fully 
turbulent. 

Table  3.  A  Summary  of  Experimental  Case  Studies  Used  for  Computational  Model  Validation 

(Adapted  from  reference  33.) 


Case  Name 

Wind  Velocity 
Magnitude 
(rms'1) 

Turbulence 

Intensity 

(%) 

Reynolds  Number 
(Re) 

R^Cylinder 

Repooi 

Low 

1.9 

10.5 

4.76  105 

2.45  106 

Medium 

5.4 

22.2 

1.35  106 

C\ 

o 

os 

High 

10.2 

16.7 

2.56  106 

1.32  107 

The  addition  of  turbulence  equations  was  an  unfortunate  but  necessary  requirement  because  not 
all  of  the  spatial  and  temporal  scales  associated  with  the  high  Re  hydrocarbon  pool  fire  flow 
environment  can  be  explicitly  modeled.  The  time-averaged  Reynolds-Averaged  Navier-Stokes 
(RANS)  equations  simulated  turbulent  flow  conditions  via  the  realizable  k-  model  [47],  This 
method  proved  most  effective  in  capturing  the  experimentally  observed  near-field  mean  flow 
structures  and  overall  heat  transfer  to  immersed  structures  while  still  maintaining  robust 
numerical  stability.  Turbulence  equations  require  auxiliary  turbulent  boundary  conditions. 
Proper  boundary  conditions  were  often  not  obvious,  and  generally  have  a  strong  effect  on 
computational  stability.  Turbulent  kinetic  energy  ( k )  and  eddy  dissipation  rate  (  )  boundary 
conditions  were  derived  from  the  TI  of  the  flow  and  large-scale  eddy  (characteristic  turbulent) 
length.  The  TI  was  extracted  from  table  3  wind  velocities  for  each  case.  The  turbulent 
characteristic  length  was  more  ambiguous,  frequently  generalized  by  a  relationship  with  the 
characteristic  length  and  the  overall  flow  domain.  An  expression  correlating  the  local  boundary 
layer  thickness  to  the  large-scale  eddy  length  produces  a  value  of  0.256  m,  which  was  used  for 
each  validation  case  [48], 

3.3.2  Combustion. 

Chemical  reaction  kinetics  were  simulated  using  the  non-premixed  combustion  model  in  which 
heated  jet  fuel  vapor  was  injected  at  the  location  of  the  horizontal  liquid  pool  surface  normal  to 
the  oncoming  turbulent  atmospheric  cross  flow.  As  the  fuel  vapor  mixed  with  the  entraining 
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oxygen,  a  fast  chemistry  or  “mixed-is-bumt”  approximation  was  used  to  determine  the  resulting 
combustion  plume  characteristics.  Also  referred  to  as  the  mixture  fraction  approach,  this  method 
employed  a  probability  density  function  to  statistically  determine  the  local  combustion  chemistry 
and  thermodynamics  based  upon  local  fuel-air  turbulent  mixing.  This  approach  was  valid  for 
most  of  the  flow  field  since  chemical  reaction  time  scales  are  much  faster  compared  to  turbulent 
fluid  motion  time  scales  present  in  pool  fire  combustion  plumes.  These  assumptions  can  lose 
validity  as  crosswinds  amplify  formulating  quasi-premixed  conditions  before  combustion 
initiates.  Figure  5  relates  the  chemical  and  physical  flow  time  scales  graphically.  Because  the 
large  spatial  scales  of  a  full-scale  aircraft  fire  environment  demand  significant  computing 
resources,  the  non-premixed  combustion  was  further  advantageous  requiring  the  solution  of  only 
two  more  transport  equations:  the  mixture  fraction  mean  and  variance  [35]. 


Fast  time  scales:  Intermediate  Slow  time  scales: 

Formation  of  time  scales:  Soot  formatoin 

radicals  CO  oxidation 


transport,  and  turbulence 

Figure  5.  Comparison  of  Physical  and  Chemical  Combustion  Time-Scales 

(Adapted  from  reference  35.) 

Energy  terms  were  treated  non-adiabatically  to  allow  for  heat  transfer  to  the  surrounding  flow 
field  and  engulfed  structures.  Inlet  diffusion  was  activated  to  allow  for  enhanced  fuel  air  mixing 
near  the  pool  surface  as  evidenced  by  the  experimental  observations.  Because  the  standard 
FLUENT  Jet  A  chemical  mechanism  was  chemically  analogous  to  JP-8  for  the  purpose  of  this 
study,  it  was  implemented  for  the  fuel  [48]. 

3.3.3  Radiation  Heat  Transfer. 

For  significantly  large-diameter  hydrocarbon  pool  fires,  radiation  heat  transfer  plays  a  critical 
role  in  overall  heat  transfer  due  to  the  emissive  power  of  the  flame  zone  participating  with  the 
local  surroundings.  Heavy  hydrocarbon  mixtures,  such  as  JP-8  and  other  diesel  derivatives, 
typically  do  not  bum  stoichiometrically  in  ambient  air  due  to  naturally  occurring  oxygen-rich 
and  oxygen-starved  currents  induced  by  the  thermally  buoyant  plume.  This  causes  inefficient 
combustion  and  in  turn  produces  heavy  soot,  which  absorbs  a  large  portion  of  the  emitted 
radiation,  as  evidenced  in  figure  2(a).  Radiation  heat  transfer  was  described  by  the  discrete 
ordinates  (DO)  model.  The  DO  model  solved  the  radiation  heat  transfer  equation  in  a  transport 
(non-ray  tracing)  format  analogous  to  the  energy  and  fluid  flow  equations  of  motion.  Radiation 
intensity  was  solved  for  using  a  preset  number  of  discrete  solid  angles  projected  through  the  flow 
domain  from  the  radiant  source  through  the  participating  medium — in  this  case,  the  air  and 
combustion-related  species  and  products.  Gray  gas  model  assumptions  were  applied  to 
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approximate  jet  fuel  fire  conditions  as  radiant  intensity  does  not  vary  significantly  with 
wavelength  [48  and  49], 

3.3.4  Soot  Formation. 

Soot  formation  plays  an  integral  role  in  determining  the  amount  of  local  energy  absorbed  during 
the  combustion  process  and  in  turn  affecting  the  total  heat  transfer  to  objects  engulfed  in  the 
flame  zone.  Soot  generation  was  governed  by  the  Moss-Brookes  method,  which  solves  soot 
nucleation  and  growth  as  a  general  transport  function  analogous  to  the  conservation  equations. 
The  Moss-Brookes  method  provided  the  most  accurate  heat  absorption  balance  to  generation 
sources  and  was  the  best  method  to  match  experimental  validation  flame  temperatures,  ft  also 
provided  a  more  qualitatively  accurate  particulate  cloud  compared  to  other  approaches  [50]. 

3.3.5  Physical  Flow  Domain. 

The  physical  flow  domain  was  designed  to  reflect  the  experimental  validation  conditions 
developed  by  SNL  and  NAWCWPNS.  To  minimize  cell  count,  a  multiblock  hybrid  mesh 
topology  was  generated  using  Pointwise  Gridgen  .  To  accommodate  multiple  crosswind  flow 
directions  and  minimize  grid  regeneration,  a  cylindrical  far-field  domain  was  used  that  allowed 
for  rotational  symmetry  about  a  central  vertical  axis.  A  hexahedral,  structured  far  field  with  high 
aspect  ratio  cells  near  the  ground  was  used  to  resolve  turbulent  boundary  layer  effects.  A 
structured  topology  also  was  used  for  the  fuel  pool  and  mock  fuselage  near  field  boundary  zones. 
An  unstructured  block  made  of  isotropic  tetrahedrons  and  pyramids  was  used  to  link  near  field 
and  far-field  structured  zones.  Approximately  1.7  million  cells  were  used  for  each  case.  Figure 
6  shows  a  near- field  schematic  with  relevant  dimensions  reflecting  experimental  schematics,  and 
shows  the  near-field  surface  mesh  of  the  mock  fuselage,  fuel  pool,  and  local  surrounding  area. 
Figure  7  shows  the  geometric  mesh  in  its  entirety  depicting  the  far-field,  outer-cylindrical 
domain  along  with  multiple  interior  cuts  describing  local  grid  topology. 
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(a)  (b) 

Figure  6.  Experimental  Validation  Case  Showing  (a)  a  Near-Field  View  of  the  Relative 
Dimensions  of  the  Fuel  Pool  and  Mock  Fuselage  and  (b)  an  Isometric  View  of  the  Equivalent 

Geometric  Surface  Mesh 


Mock  Fuselage  and  Ground 
Grid  Cross-Section 


Mockup  Fuselage  and 
Fuel  Pool  Near  Field 


Figure  7.  Experimental  Validation  Case  Geometric  Mesh 
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3.3.6  Boundary  Conditions. 


Boundary  and  operating  conditions  reflected  the  measurements  recorded  during  the  experimental 
validation  where  applicable,  with  some  adjustments  made  for  simplicity.  To  define  the 
prevailing  wind  speed  and  direction,  the  cylindrical  far-field  boundary  was  divided  into  three 
surfaces.  The  windward  far-field  boundary  was  defined  as  a  constant  velocity  inlet  using  the 
mean  measured  wind  speed.  The  opposing  leeward  far-field  boundary  and  the  ceiling  were 
defined  as  a  pressure  outlet  to  balance  the  flow.  The  fuel  pool  surface  was  treated  as  a  gaseous 
mass  flow  inlet  set  at  the  mean  temperature  of  508  K  recorded  by  experiments  [33],  Fuel  vapor 
was  injected  normal  to  the  domain  at  16.46  kg^'1,  a  value  extrapolated  from  the  average  fuel 
recession  rate  and  pool  geometry.  The  mock  fuselage  and  surrounding  ground  area  were 
modeled  as  infinitely  thin,  nonadiabatic  walls.  All  bounded  surfaces  besides  the  fuel  pool  fire 
inlet  were  initialized  to  300  K  to  simulate  start-up  conditions.  As  steady-state  combustion 
conditions  equilibrated  during  the  solution  process,  convective  and  radiation  heat  transfer  forced 
the  mock  fuselage,  the  flame  zone,  and  the  surrounding  ground  area  to  heat  up.  Fuel  pool  mass 
flow  inlet  conditions  also  floated  to  accommodate  feedback  from  the  flow  domain,  although 
steady-state  conditions  stayed  very  close  to  the  prescribed  conditions.  Figure  8  shows  the 
location  and  boundary  conditions  enforced  relative  to  the  far  field.  Figure  9  depicts  boundary 
conditions  from  the  near  field  perspective. 


Prevailing 
Wind  Direction 


Velocity  Inlet 

Pressure  Outlet 

Mockup  Fuselage  Wall 


"J 


Mass  Flow  Inlet 


Y  (Out  of  Paper) 


Figure  8.  Far-Field  Schematic  Outlining  CFD  Boundary  Types  and  Conditions 
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Pressure  Outlet  Atmospheric  Cross  Wind  Velocity  Inlet 

Pre-Set  @  300  K  Wind  Speed  /  Direction  Dependent  Pre-Set  @  300  K 


Pre-Set  @  300  K  Pre-Set  @  300  K 


Figure  9.  Near-Field  Schematic  Outlining  CFD  Boundary  Types  and  Conditions 
3.3.7  Numerics. 

Model  simulations  were  executed  using  an  implicit,  segregated  pressure -based  solver  using  the 
SIMPLE  (Semi-Implicit  Method  for  Pressure-Linked  Equations)  algorithm  for  pressure-velocity 
coupling  [48].  Steady-state  simulations  were  solved  with  second-order  spatial  accuracy,  with 
time-accurate  solutions  requiring  first-order  accuracy  enabled  for  radiation  transport  to  maintain 
solution  stability.  Temporal  discretization  was  maintained  using  the  bounded  second-order 
implicit  method.  Time-accurate  calculations  were  solved  with  a  fixed  time  step.  Although 
absolute  convergence  criteria  were  not  enforced,  residuals  were  observed  to  be  on  the  order  of 
10"4  or  lower. 

3.4  RESULTS. 

Computational  validation  cases  were  executed  on  an  8-processor  Harpertown  Intel®  Xeon® 
workstation  taking  approximately  50  central  processing  unit  (CPU)  hours  per  case,  or  about  6 
hours  real  time.  All  reproduced  experimental  data  used  for  comparison  was  time-averaged  over 
a  finite  time  period  as  outlined  in  table  1.  CFD  model  validation  studies  were  solved  steady- 
state.  Mock  fuselage  surface  pressure,  heat  flux,  temperature,  and  flame  temperature  were  the 
data  sets  compared. 

3.4.1  Surface  Pressures. 

DPGs  were  installed  on  the  mock  fuselage  surface  to  record  the  change  in  pressure  between  the 
windward  stagnation  surface  and  leeward  trailing  surface.  Measurements  were  only  available  for 
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the  medium  wind  condition  case  [33].  Figure  10  shows  the  unsteady  experimental  (EXP)  DPG 
measurements  for  the  mock  fuselage  (MF)  west,  center,  and  east  locations,  as  well  as  the  relative 
fuel  pool  (Flame)  surface  to  far- field  DPG  described  in  section  3.2.  The  specific  quasi-steady 
time  period  considered  for  data  was  recorded  between  the  6  min  40  s  and  9  min  36  s  mark. 
Figure  1 1  shows  the  equivalent  CFD  static  pressure  calculations  for  the  west,  center,  and  east 
locations  reported  at  each  discrete  computational  cell  surface  location  around  the  mock  fuselage 
circumference.  The  windward  and  leeward  locations  in  addition  to  the  fuselage  top  and  bottom 
are  highlighted  with  a  black  circle.  Their  locations  relative  to  the  mock  fuselage  circular  cross- 
section  is  shown  in  the  figure’s  lower-left  comer.  An  average  CFD  differential  pressure  of  about 
55  Pa  was  reported  for  all  3  fuselage  surface  locations  showing  good  agreement  with  the 
observed  results.  The  CFD  model  calculated  an  approximate  12  Pa  differential  pressure  for  the 
flame,  which  also  compared  well  to  the  experiment  from  figure  10.  Figures  12  and  13  show  the 
CFD  total  pressure  and  velocity  vectors,  respectively,  in  the  mock  fuselage  near  field  along  a 
centerline  cross-section.  The  mock  fuselage  wake  was  asymmetric  due  to  a  combination  of  its 
proximity  to  the  ground  coupled  with  the  thermal  buoyancy  generated  near  the  fuel  surface. 
Velocity  vectors  show  a  similar  trend  with  corresponding  local  flow  directions  depicting  the 
wake  region  characteristic  of  bluff-body  shedding. 


CFD  =  Computational  Fluid 
Dynamic 

EXP  =  Experimental 
MF  =  Mock  Fuselage 


Figure  10.  Medium  Wind  Condition  Transient  Experimental  DPG  Measurements  on  the  Mock 
Fuselage  Compared  to  Time-Averaged  Steady  CFD  Results  (Adapted  from  reference  33.) 
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M  ock  Fuselage  Axial  Cross-Section  Location  (m) 


Figure  11.  The  CFD  Static  Pressure  Calculated  About  the  Mock  Fuselage  Cylinder 
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Contours  of  Total  Pressure  (pascal) 


Figure  12.  Centerline  CFD  Mock  Fuselage  Near-Field  Total  Pressure  Contours  of  the  Medium 

Wind  Condition 
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Verocity  Vectors  Coioied  By  Velocity  Magnitude  (m/s) 


Figure  13.  The  CFD  Mock  Fuselage  Near-Field  Velocity  Magnitude  Vectors  of  the  Medium 

Wind  Condition 


3.4.2  Flame  Temperatures. 

Experimental  flame  temperature  measurements  from  the  case  studies  outlined  in  table  3  were 
recorded  and  time-averaged  using  contour  plots  in  three  equally  spaced  planes  perpendicular  to 
the  fuel  pool  in  figures  14,  15,  and  16.  Experimental  validation  plots  are  aligned  on  the  left,  and 
the  corresponding  CFD  plots  are  aligned  on  the  right.  Figure  legends  denoting  each  data  plane 
are  located  to  the  far  left  on  each  plot.  (See  the  dashed  outlined  planes  from  figure  3(b)  for  the 
detailed  plane  reference.)  Plane  measurements  were  approximated  because  all  the  TCs  did  not 
lie  in  the  same  exact  plane.  The  intersection  with  the  mock  fuselage  is  represented  by  the  black 
circular  cross-section  in  each  plot.  Measured  and  calculated  temperatures  ranged  from  300  K 
(ambient)  to  upwards  of  1600  K  within  and  near  the  flame  zone.  CFD  color  maps  were 
converted  to  gray  scale  and  discretized  to  reflect  the  same  ranges  reported  experimentally  [33], 

Computational  results  compared  well  capturing  similar  magnitudes  and  profile  trends.  Similar 
agreement  was  illustrated  for  the  low  crosswind  case,  with  the  high  crosswind  case  showing  the 
greatest  divergence  most  prevalent  along  the  leeward  side  of  the  fuselage  and  adjacent  trailing 
wake  zone.  Heat  transfer  was  under  predicted  significantly  in  this  regime.  Because  fuselage 
surface  differential  pressure  data  was  unavailable  for  the  high  crosswind  condition,  it  was 
difficult  to  speculate  how  the  model  broke  down.  From  inspection,  it  was  theorized  turbulent 
flow  separation  was  predicted  prematurely  creating  a  lifted  flame  condition  and  in  turn 
minimizing  the  heat  transfer  to  the  leeward  surface. 

Figure  14  shows  good  agreement  between  experimental  and  CFD  results  for  the  low  crosswind 
case.  The  low  crosswind  case  maintained  the  most  acute  angle  with  respect  to  the  mock  fuselage 
position  and  depicted  the  greatest  asymmetry  in  flow  between  each  plane.  The  low  crosswind 
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case  also  provided  the  lowest  TI,  or  lowest  ratio  of  wind  fluctuations  relative  to  the  mean  flow 
creating  the  calmest  experimental  wind  condition.  Still  air  to  low  crosswind  conditions  create 
the  best  conditions  for  a  diffusion-dominated  flame,  reflecting  an  environment  most  suitable  for 
the  mixture  fraction  combustion  approach  taken.  Figure  15  shows  similar  agreement  for  the 
medium  wind  case.  The  medium  wind  case  had  the  most  relative  normal  angle  with  respect  to 
the  mock  fuselage  position,  which  created  the  strongest  symmetry  among  all  cross-sections. 
This  is  visibly  evident  from  the  plots.  Medium  crosswind  conditions  generated  the  largest 
crosswind  TIs  of  all  three  cases  and  provided  the  strongest  flow  mixing,  which  created  the  most 
favorable  model  conditions  for  RANS  calculations.  Figure  16  shows  comparisons  between  the 
high  crosswind  cases  showing  significant  disparity  in  leeward  wake  region  flame  temperature 
magnitudes.  Here,  two  theories  are  postulated  based  upon  employed  model  assumptions. 
Convective  forces  due  to  enhanced  crosswind  conditions  created  greater  flow  separation  around 
the  leeward  side  of  the  cylinder.  With  the  employment  of  a  simplified  RANS  turbulence  model, 
flow  separation,  detachment,  and  downstream  recirculation  became  increasingly  hard  to  predict. 
Since  recirculation  enhances  fuel-air  mixing  and  elevates  heat  transfer,  poor  flow  prediction  in 
this  area  could  affect  thermal  predictions.  Further,  higher  crosswind  forces  challenge  the 
assumptions  made  by  the  mixture  fraction  combustion  modeling  approach  conducive  to  diffusion 
flames.  With  an  increase  in  forced  air  entrainment,  heated  fuel  vapors  have  a  longer  residence 
time  to  premix  with  a  stronger  prevailing  wind  before  combusting.  This  can  delay  ignition 
sweeping  high  heat  release  zones  away  from  highly  turbulent  fuel-air  mixing  zones. 
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Figure  14.  Low  Crosswind  Case  Showing  Flame  Temperature  Contour  Comparisons 
(Experimental  data  was  adapted  from  reference  33.) 
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Figure  15.  Medium  Crosswind  Case  Showing  Flame  Temperature  Contour  Comparisons 
(Experimental  data  was  adapted  from  reference  33.) 
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Figure  16.  High  Crosswind  Case  Showing  Flame  Temperature  Contour  Comparisons 
(Experimental  data  was  adapted  from  reference  33.) 

3.4.3  Surface  Temperatures. 

Temperatures  were  recorded  along  the  mock  fuselage  surface  similar  to  the  way  in  which  flame 
temperature  measurements  were  taken  in  section  3.4.2,  but  instead  using  surface  mount  TCs  as 
opposed  to  standing  rakes.  Experimental  measurements  are  again  aligned  on  the  left,  and  CFD 
results  are  aligned  on  the  right.  Mock  fuselage  surfaces  are  broken  down  to  the  windward  and 
leeward  view,  annotated  by  the  legend  on  the  far  left.  Experimental  circular  fuselage  surface 
data  is  projected  onto  a  flat  surface  defined  via  degrees,  and  CFD  skin  data  is  presented  using  an 
orthonormal  perspective.  A  similar  scale  from  300  K  to  about  1550  K  was  implemented  to 
depict  surface  temperatures. 

Figures  17,  18,  and  19  report  mock  fuselage  surface  temperatures  for  the  low,  medium,  and  high 
wind  speed  cases,  respectively.  Local  surface  temperatures  for  all  three  cases  reflect  similar 
results  reported  in  the  mock  fuselage  near  field  local  flame  temperatures  documented  in 
section  3.4.2.  Surface  temperatures  show  good  agreement  for  the  low  and  medium  crosswind 
cases,  with  significant  divergence  in  comparing  the  high  wind  speed  experimental  and  CFD  case 
for  the  same  postulated  reasons  mentioned  at  the  end  of  section  3.4.2. 
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Figure  17.  Low  Crosswind  Case  Showing  Mock  Fuselage  Surface  Temperature  Contours 
Comparisons  (Experimental  data  was  adapted  from  reference  33.) 
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Figure  18.  Medium  Crosswind  Case  Showing  Mock  Fuselage  Surface  Temperature  Contour 
Comparisons  (Experimental  data  was  adapted  from  reference  33.) 
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Figure  19.  High  Crosswind  Case  Showing  Mock  Fuselage  Surface  Temperature  Contour 
Comparisons  (Experimental  data  was  adapted  from  reference  33.) 

3.4.4  Surface  Heat  Fluxes. 

Heat  flux  was  recorded  along  the  mock  fuselage  surface  similar  to  surface  temperatures 
mentioned  in  section  3.4.3  and  appears  in  figures  20,  21,  and  22.  The  same  data  presentation 
methods  were  used.  A  single  scale  from  approximately  5  to  230  kWxrf  was  used  for  the  low 
crosswind  case  in  figure  20,  and  40  to  310  kW>m"2  was  used  for  the  medium  and  high  crosswind 
cases  to  clarify  the  higher  recorded  heat  fluxes  at  higher  crosswind  conditions  in  figures  21  and 
22,  respectively.  Surface  heat  flux  showed  significant  gains  experimentally  as  crosswind 
increased  for  the  low  and  medium  wind  speed  cases  and  diverged  for  the  high  wind  speed  case. 
Similar  conclusions  were  reached  for  heat  flux  comparisons  as  the  flame  and  surface 
temperatures. 


24 


Figure  20.  Low  Crosswind  Case  Showing  Mock  Fuselage  Surface  Heat  Flux  Contour 
Comparisons  (Experimental  data  was  adapted  from  reference  33.) 
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Figure  21.  Medium  Crosswind  Case  Showing  Mock  Fuselage  Surface  Heat  Flux  Contour 
Comparisons  (Experimental  data  was  adapted  from  reference  33.) 
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Figure  22.  High  Crosswind  Case  Showing  Mock  Fuselage  Surface  Heat  Flux  Contour 
Comparisons  (Experimental  data  was  adapted  from  reference  33.) 

3.4.5  Flow  Visualization. 

Color  contour  temperature  plots  showing  the  fuel  pool  surface,  mock  fuselage,  and  three  cross- 
sectional  planes  are  depicted  for  all  three  crosswind  cases  in  figure  23.  Low  crosswinds  exhibit 
the  lowest  temperatures  due  to  minimal  fuel-air  mixing  combined  with  the  most  acute  wind 
direction  with  respect  to  the  mock  fuselage.  Medium  crosswinds  had  the  largest  wind 
fluctuations  supporting  greater  fuel-air  mixing.  High  crosswinds  had  the  least  amount  of  wind 
fluctuations  but  the  strongest  mean  flow  creating  a  more  aggressively  stretched  wake  region. 
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Figure  23.  Flame  and  Surface  Contour  Temperatures  for  the  Case  CFD  Validation  Case  Studies 
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4,  AIRCRAFT  COMPUTATIONAL  ANALYSIS. 


4.1  OVERVIEW. 

Three  aircraft  geometries  were  selected  to  analyze  relative  thermal  responses  to  similar  jet  fuel 
pool  fire  conditions.  Low  and  moderate  crosswind  conditions  with  10%  free-stream  TI  were 
imposed  in  the  direction  normal  to  each  aircraft  fuselage  to  approximate  atmospheric  conditions. 
These  settings  were  selected  based  upon  good  agreement  with  experimental  data  for  those 
respective  flow  regimes.  The  B-707-320B  was  selected  to  represent  the  physical  features  of  a 
traditional  airliner  in  terms  of  size  and  shape.  The  B-777-200ER  was  one  of  two  NLA  designs 
chosen  due  to  its  large  frame  classification  in  addition  to  its  circular  cross-sectional  fuselage  and 
twin-engine  build.  The  third  airframe  analyzed  was  the  NLA  A380-800  equipped  with  four 
engines  and  uniquely  identified  by  its  elliptically  shaped  hull.  For  brevity,  aircraft  will  be 
referred  to  by  their  family  name. 

The  National  Fire  Protection  Association  (NFPA)  403  practical  critical  area  (PCA)  fire  scenario 
was  the  standard  chosen  for  comparison  because  it  provided  a  well-established  geometric  basis 
to  generate  similar  scale  aircraft  fires.  It  is  formulated  from  the  following  relationship: 

PCA  =  0.67  X  [L  X  (K  +  W)]  (1) 

where  L,  W,  and  K  are  the  aircraft  length,  fuselage  width,  and  a  constant  based  upon  L, 
respectively  [1].  Figure  24  relates  aircraft  position  to  PCA.  Table  4  provides  relative  area 
statistics  of  each  aircraft  fire  environment,  and  table  5  summarizes  aircraft  fuel  capacity  and 
passenger  limitations  dependent  upon  internal  cabin  configuration  [51-53].  It  should  be  noted 
here  that  aircraft  PCA,  fuel,  and  passenger  capacity  do  not  scale  proportionally.  Figures  25  and 
26  show  multiple  views  of  each  aircraft  and  their  unique  PCA  scaled  relative  to  one  another. 
The  same  geometries  were  also  used  for  CFD  physical  domain  construction. 


Figure  24.  The  PCA  Size  and  Location  With  Respect  to  the  Aircraft 
(Adapted  from  reference  1 .) 
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Table  4.  A  Summary  of  Aircraft  Wetted  Surface  and  PCA  Ratios 


Aircraft 

Surface  Area 
(m2) 

PCA 

(m2) 

Ratio 

(-) 

B-707 

1252 

1054 

1.19:1 

B-777 

2341 

1545 

1.51:1 

A3  80 

3975 

1784 

2.23:1 

Table  5.  A  Summary  of  Aircraft  Fuel  Capacity  and  Passenger  Limitations  [51-53] 


Aircraft 

Fuel  Capacity 
(1) 

Passenger  Limit* 

B-70715" 

90,290 

147  or  202 

B-77715" 

171,170 

301  or  400 -440 

A380L55J 

323,546 

525  or  853 

*Configuration  Dependent 


Top  View 


Figure  25.  A  Comparison  of  the  Relative  Scale  and  Shape  Between  Each  Aircraft 
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Isometric  View 

Figure  26.  An  Alternative  Perspective  Comparing  the  Relative  Scale  and  Shape  Between 

Each  Aircraft 


4.2  PHYSICAL  DOMAIN. 

Flow  domains  for  all  three  aircraft  were  composed  similarly  to  the  computational  validation 
studies.  A  maximum  average  surface  mesh  spacing  of  10  cm  was  applied  to  each  aircraft,  with 
areas  of  high  curvature  more  finely  resolved.  Mesh  sizes  ranged  from  9.5  (B-707)  to  10.9 
(A380)  million  cells.  Figure  27  depicts  various  aspects  of  the  A380  geometric  mesh. 
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Figure  27.  The  A380  Geometric  Mesh 

4.3  BOUNDARY  CONDITIONS. 

Boundary  conditions  for  each  aircraft  analysis  were  imposed  similarly  to  the  validation  studies. 
The  aircraft  surface  was  treated  in  the  same  manner  as  the  mock  fuselage  surface,  and  the  PCA 
was  treated  in  a  similar  manner  as  the  pool  fire  zone  by  extrapolating  the  fuel  recession  mass 
flow  rate  to  compensate  for  the  increase  in  area.  Free-stream  conditions  and  surrounding  ground 
area  were  treated  identically. 

4.4  RESULTS. 

Steady-state  as  well  as  time-accurate  simulations  using  a  50-ms  fixed  time  step  over  a  120-s 
period  were  calculated  for  each  of  the  three  aircraft  exposed  to  low  (2  rms'1)  and  medium 
(5  rms'1)  crosswind  conditions.  Unsteady  calculations  were  initialized  from  the  steady-state 
solution  allowing  30  s  of  transition  to  fully  developed  unsteady  flow.  All  cases  were  executed 
on  a  Department  of  Defense  Supercomputing  Resource  Center  SGI®  Altix®  ICE.  Hardware 
utilization  consisted  of  64  Nehalem  Intel®  Xeon®  processor  cores  connected  via  InfiniBand™. 
Each  simulation  took  approximately  4000  CPU  hours  (about  60  real-time  hours)  to  complete. 
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4.4,1  Surface  Temperatures  and  Heat  Fluxes. 


Aircraft  surface  thermal  response  was  a  critical  factor  recorded  for  each  time  step  to  interpret 
dynamic  aircraft-fire  interaction.  Figure  28  shows  transient  mean  temperatures  and  heat  fluxes 
averaged  over  the  entire  aircraft  surface  during  a  fully  developed  60-s  unsteady  period  exposed 
to  low  and  medium  crosswinds.  Results  indicated  that  low  crosswinds  allow  the  periodic  mean 
fire  plume  pulsations  to  dominate  heat  transfer  to  the  aircraft.  Conversely,  medium  crosswind 
conditions  force  a  higher  rate  of  aircraft-plume  interaction  resulting  in  faster  secondary  flow 
structures.  This  causes  increased  turbulent  fuel-air  mixing  and  greater  heat  transfer.  Table  6 
shows  a  quantitative  comparison  of  the  mean  thermal  response  between  each  aircraft;  this 
comparison  shows  heat  transfer  magnitudes  approximately  doubling  between  low  and  medium 
crosswinds  conditions,  a  trend  established  during  the  CFD  validation  study.  Percent  differences 
are  calculated  relative  to  the  B-707.  It  is  further  speculated  that  as  crosswind  magnitudes 
increase,  overall  heat  transfer  becomes  increasingly  independent  of  secondary  changes  in  aircraft 
geometry  (fuselage  shape,  engine  layout,  wing  planform,  etc.)  as  long  as  the  mean  wake  structure 
created  by  the  aircraft  body  is  not  significantly  altered. 


B707  Temperature 

B707  Heat  Flux 

-  B777  Temperature 

.  B777  Heat  Flux 

-  A380  Temperature 

. A380  Heat  Flux 

Figure  28.  Comparisons  of  Mean  Aircraft  Surface  CFD  Temperatures  and  Heat  Fluxes  Exposed 

to  (a)  Low  and  (b)  Medium  Crosswinds 
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Table  6.  A  Summary  of  Time- Averaged  Mean  Aircraft  Surface  Thermal  Signatures 


Aircraft 

Mean  ±  RMS 
Temperature 
(K) 

Difference 

(%) 

Mean  ±  RMS 
Heat  Flux 
(kWxn2) 

Difference 

(%) 

Low  Crosswind 

B-707 

613  ±25 

— 

16.8  ±2.9 

— 

B-777 

701  ±42 

14.4 

31.4  ±7.5 

86.9 

A3  80 

639  ±  26 

4.3 

22.4  ±3.3 

33.3 

Medium  Crosswind 

B-707 

1050  ±65 

— 

59.7  ±4.9 

— 

B-777 

957  ±  82 

8.6 

55.6  ±6.5 

6.9 

A3  80 

949  ±  58 

9.6 

52.3  ±5.1 

12.4 

Figures  29  through  32  illustrate  time-averaged  aircraft  surface  temperatures  and  heat  fluxes 
exposed  to  medium  crosswind  conditions  for  a  60-s  period.  Aircraft  fuselage  lengths  are 
normalized  for  easier  interpretation  of  results.  Mean  temperatures  peaked  at  approximately 
1200  K,  and  mean  heat  fluxes  peaked  at  about  100  kWxn'2.  RMS  temperatures  fluctuated 
between  10%  and  70%  of  the  time-averaged  value.  Low  crosswind  results  depicted  similar 
contour  profile  trends  but  at  lower  magnitudes,  as  quantified  in  table  6. 
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Figure  29.  Top  View  of  Time- Averaged  Aircraft  Surface  Temperatures  Exposed  to 

Medium  Crosswinds 
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Figure  30.  Bottom  View  of  Time- Averaged  Aircraft  Surface  Temperatures  Exposed  to 

Medium  Crosswinds 
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Figure  3 1 .  Top  View  of  Time- Averaged  Aircraft  Surface  Heat  Flux  Exposed  to 

Medium  Crosswinds 
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Contours  of  Surface  Heat  Flux  (W/m2) 


Figure  32.  Bottom  View  of  Time- Averaged  Aircraft  Surface  Heat  Flux  Exposed  to 

Medium  Crosswinds 


Certain  locations  exhibited  strong  three-dimensional  flow  effects  such  as  around  engine  nacelles, 
the  fuselage  nose,  and  horizontal  and  vertical  stabilizers.  These  locations  that  also  existed  in  the 
flame  zone  generated  enhanced  heat  transfer  due  to  greater  turbulent  fuel-air  mixing  and  more 
exposure  to  entrained  air.  Areas  in  close  proximity  to  the  ground  and  central  to  the  body  were 
exposed  primarily  to  more  oxygen-starved,  fuel-rich  conditions  reducing  local  heat  transfer. 

Near-field  time-averaged  aircraft  surface  and  PCA  fuel  pool  temperatures  are  shown  in  figures 
33  and  34,  which  depict  the  time-averaged  effect  of  the  aircraft-fire  plume  interaction  on  the 
surrounding  area.  Enhanced  areas  of  heat  transfer  along  fuel-air  interfaces  are  identified.  Three- 
dimensional  heat  transfer  effects  are  also  more  easily  observed.  Increases  in  surface  temperature 
in  the  near-field  ground  area  are  attributed  to  the  crosswind  flow  constrictions  induced  by  aircraft 
proximity  with  the  ground.  Because  the  B-707  maintains  the  lowest  overall  average  ground 
clearance  of  all  three  aircraft,  local  fuel-air  mixing  in  these  regions  is  enhanced  causing  the 
largest  rise  in  temperature. 
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Figure  33.  Time- Averaged  Near-Field  Aircraft  Surface  Fuel  Zone  and  Ground  Temperatures 
Exposed  to  Medium  Crosswinds  (Isometric  View) 
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Figure  34.  Time-Averaged  Near-Field  Ground  and  Fuel  Surface  Temperatures  Exposed  to 

Medium  Crosswinds  (Top  View) 


4.4.2  Flame  Temperatures. 

Figures  35  and  36  show  time-averaged  flame  temperatures  from  both  far-  and  near-field 
perspectives,  respectively.  Radial  temperature  profiles  emanating  from  the  fire  plume  center 
compared  well  among  all  three  aircraft,  increasing  in  width  as  PCA  increased  and  showing  peak 
flame  temperatures  of  approximately  1000  K.  Figure  35  also  shows  mean  flame  height  also 
increased  with  the  base  fire  size. 
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Figure  35.  Time-Averaged  Far-Field  Fire  Plume  Temperatures  Exposed  to  Medium  Crosswinds 

(Front  View) 
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Boeing  707 
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Figure  36.  Time-Averaged  Near-Field  Fire  Plume  Temperatures  Exposed  to  Medium 

Crosswinds  (Front  View) 


4.4.3  Flame  Velocities. 

Figure  37  shows  the  same  reference  frame  as  figure  34  illustrating  far-field  fire  plume  velocity 
magnitudes  for  each  aircraft  exposed  to  medium  crosswinds.  Similar  to  flame  height  and  width, 
core  flame  velocities  increased  with  increasing  PCA,  and  peak  magnitudes  of  approximately 
35  nDS1  were  calculated. 
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Figure  37.  Time-Averaged  Far-Field  Fire  Plume  Velocity  Magnitudes  Exposed  to  Medium 

Crosswinds  (Front  View) 


5.  CONCLUSIONS. 

A  computational  aircraft  pool  fire  model  was  developed  capable  of  predicting  the  major 
combustion  plume  features  and  heat  transfer  magnitudes  to  immersed  objects  exposed  to  low  and 
medium  atmospheric  crosswind  conditions.  Computational  fluid  dynamic  (CFD)  results 
supported  experimental  measurements  that  as  wind  velocity  increases,  turbulent  fuel-air  mixing 
is  enhanced,  which  results  in  increased  heat  transfer  to  bodies  engulfed  in  the  flame  zone. 
Divergence  from  experimental  data  at  high  crosswinds  can  likely  be  attributed  to  multiple  factors 
involving  the  breakdown  of  turbulence  and  combustion  model  assumptions.  Discrepancies  at 
lower  crosswind  conditions  are  likely  attributed  to  inconsistencies  related  to  quasi-steady  state 
flow  assumptions  along  with  simplifications  made  in  generating  experimental  flame  zone  planes 
for  comparison  to  CFD.  To  increase  accuracy,  higher-fidelity  Large-Eddy  Simulation  (LES) 
turbulence  models  were  examined,  but  they  exhibited  little  improvement  over  mean  and  root 
mean  square  thermal  signatures  already  calculated  and  well  validated  using  Reynolds- Averaged 
Navier-Stokes  (RANS).  Numerical  stability  also  decreased  for  LES  calculations,  decreasing 
overall  robustness  to  the  solution  process. 
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Time-averaged  aircraft  analysis  suggested  a  consistent  correlation  between  different  aircraft 
geometries  exposed  to  the  same  crosswind  conditions  when  considering  fire  scenarios 
proportional  to  their  overall  size.  As  crosswind  increased,  aircraft  thermal  response  became  less 
dependent  upon  aircraft  geometry  due  to  inertial  atmospheric  forces  driving  turbulent  fuel-air 
mixing.  For  weaker  crosswinds,  aircraft-fire  plume  interaction  governed  heat  transfer.  Flame 
dynamics  showed  similar  trends  when  comparing  fire  plume  velocities  and  temperatures.  Other 
reasons  for  minor  thermal  response  differences  may  be  attributed  to  fire  scale  similarity 
misrepresentations  by  the  practical  critical  area  (PCA).  This  becomes  plausible  when  other 
characteristics,  such  as  the  aircraft  wetted  surface  area,  scale  more  linearly  with  factors  like  fuel 
and  passenger  capacity. 

Overall  interpretation  of  the  data  infers  nominal  design  changes  in  aircraft  geometry  (fuselage 
shape,  wing  planform,  and  engine  layout)  do  not  create  an  extraordinary  fire  protection  threat 
when  exposed  to  similar  scale  flame  conditions.  This  holds  true  for  both  New  Large  Aircraft 
(NLA)  as  well  as  traditionally  sized  aircraft.  Although  unique  instantaneous  heat  signatures 
were  observed,  general  trends  remained  the  same.  All  thermal  attributes  stayed  within  the  same 
order  of  magnitude  and,  in  the  majority  of  instances  varied  less  than  15%,  most  notably  when 
exposed  to  higher  heat  transfer  conditions.  However,  with  NLA  carrying  larger  amounts  of  fuel 
compared  to  the  majority  of  commercial  passenger  carrying  aircraft,  NLA  have  the  capacity  to 
create  a  significantly  larger  accidental  fire  threat. 

As  hybrid  RANS-LES  turbulence  modeling  techniques  gain  widespread  reliability,  they  could  be 
considered  to  improve  overall  accuracy  of  the  simulation  process.  Novel  meshing  techniques, 
like  anisotropic  boundary  layers  to  simplify  grid  topologies  along  with  other  mesh-reducing 
optimizations,  could  also  be  implemented  as  their  dependability  increases  to  reach  faster  solution 
turnaround  times.  Multiscale  fire  scenarios  could  also  be  considered  for  analysis,  as  sole  reliance 
on  PCA  similarity  could  be  suspect.  In  addition,  crosswind  magnitudes  and  directions 
representative  of  historical  fire  emergency  scenarios  could  be  introduced  to  integrate  real-world 
metrics  into  the  analysis. 
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